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We present an analysis of the thermalization rate of Higgsinos and winos based 
on the imaginary part of the two-point Green function in the unbroken phase 
of the MSSM. We use improved propagators including resummation of hard 
thermal loops and the thermalization rate is computed at the one-loop level 
in the high temperature approximation. We find that the damping is typically 
dominated by scattering with gauge bosons, resulting in a damping rate of about 
jjj ~ 0.025T, 7^ ~ 0.065T. The contribution from scattering with scalars is 
relatively small. Implications for baryogenesis are also discussed. 
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The presence of unsuppressed baryon number violating processes at high temper- 
atures within the Standard Model (SM) of weak interactions makes the generation of 
the baryon number at the electroweak scale an appealing scenario [Q. The baryon 
number violating processes also impose severe constraints on models where the baryon 
asymmetry is created at energy scales much higher than the electroweak scale [fjj. 
Unfortunately, the electroweak phase transition is too weak in the SM |5) (although 
the existence of a primordial hypermagnetic field could improve the situation ||). 
This means that the baryon asymmetry generated during the transition would sub- 
sequently be erased by unsuppressed sphaleron transitions in the broken phase. The 
most promising and well-motivated framework for electroweak baryogenesis beyond 
the SM seems to be supersymmetry (SUSY). Electroweak baryogenesis in the frame- 
work of the Minimal Supersymmetric Standard Model (MSSM) has attracted much 
attention in the past years, with particular emphasis on the strength of the phase 
transition || and the mechanism of baryon number generation || [7], || ||, [TJJ [H| . 



Recent analytical |T2|, [Uj and lattice computations |14j] have revealed that the phase 
transition can be sufficiently strongly first order if the ratio of the vacuum expectation 
values of the two neutral Higgses tan/3 is smaller than ~ 4. Moreover, taking into 
account all the experimental bounds including those coming from the requirement of 
avoiding dangerous color breaking minima, the lightest Higgs boson should be lighter 
than about 105 GeV, while the right-handed stop mass might be close to the present 
experimental bound and should be smaller than, or of the order of, the top quark mass 

Moreover, the MSSM contains additional sources of CP-violation besides the CKM 
matrix phase. However, an acceptable baryon asymmetry from the stop sector may 
only be generated through a delicate balance between the values of the different soft 
supersymmetry breaking parameters contributing to the stop mixing parameter, and 
their associated CP- violating phases ||. As a result, the contribution to the final 
baryon asymmetry from the stop sector turns out to be negligible. On the other hand, 
charginos and neutralinos may be responsible for the observed baryon asymmetry in 
the MSSM|| [D], |T5| . If the strength of the electroweak phase transition is enhanced 
by the presence of some new degrees of freedom beyond the ones contained in the 
MSSM, e.g. some extra standard model gauge singlets, light stops (predominantly the 
right-handed ones) and charginos/neutralinos are expected to give quantitatively the 
same contribution to the final baryon asymmetry. 

CP- violating sources in supersymmetric baryogenesis are more easily built up if the 
degrees of freedom in the stop and gaugino / neutralino sectors are nearly degenerate in 
mass This resonant behaviour is associated with the possibility of absorption (or 
emission) of Higgs quanta by the propagating supersymmetric particles. For momenta 
of the order of the critical temperature, this can only take place when, for instance, 



the Higgsinos and gauginos do not differ too much in mass. By using the Uncertainty 
Principle, it is easy to understand that the width of this resonance is expected to 
be proportional to the thermalization rate of the particles giving rise to the baryon 
asymmetry 0. 

From all these considerations, it is clear that the computation of the thermalization 
rate of the particles responsible for supersymmetric baryogenesis represents a neces- 
sary and crucial step towards the final computation of the baryon number. A detailed 
computation of the thermalization rate of the right-handed stop from the imaginary 
part of the two-point Green function has been recently performed in \T& by making 
use of improved propagators and including resummation of hard thermal loops. The 
thermalization rate has been computed exactly at the one-loop level in the high tem- 
perature approximation as a function of the plasma right-handed stop mass mj (T) 
and an estimate for the magnitude of the two-loop contributions which dominate the 
rate for small mj R (T) was also given. If mj (T) > T, the thermalization is dictated by 
the one-loop thermal decay rate which can be larger than T [16| . For smaller values of 
mj R (T), when the thermalization is dominated by two- loop effects (i.e. scattering), 
may be as large as 1CT 3 T [jnj. With such value, our derivative expansion is perfectly 
justified since the wall thickness can span the range (10 — 100)/T. 

The goal of this paper is to present a computation of the thermalization rate of 
Higgsinos and winos from the imaginary part of the two-point Green function in the 
unbroken phase of the MSSM. We use improved propagators including resummation 
of hard thermal loops and the thermalization rate is computed at the one-loop level 
in the high temperature approximation. 

There are two types of diagram contributing to the damping rate of charginos and 
neutralinos to one loop. The internal loop has either a scalar or a vector boson, apart 
from the fermion. We expect the leading processes to involve massless gauge bosons, 
since these diagrams are IR divergent without plasma resummation. 

Let us begin with performing an analysis for the charged H ^ . We first calculate the 
gauge boson contribution (as accurately as possible with the present understanding of 
the IR sensitivity) and then compare it with a typical diagram involving a scalar (the 
stop). We find that scattering with scalars provides typically much less damping. The 
contribution from unbroken gauge bosons is in some sense the most difficult diagram 
since it requires subtle resummation schemes. On the other hand, this is also the most 
studied case and the result can be found in the literature. The relevant diagram is 



B or W ±fi 




H±,o 



Depending on the assumptions of the mass and momentum of the external fermion 
there are a number of different cases: 

• Heavy m ^> T fermion at rest p = 0. 

In this case thermal fermions are not present but only the gauge bosons. The 



result for the damping rate 7 is [17 



Here Cf denotes the quadratic Casimir invariant of the fermion representation. 
Massless fermion at rest p — 0. 

This calculation is more involved than the heavy fermion case since one needs to 



resum all propagators and vertices. In [[18] this is done for a few different gauge 
groups with the result that 

1 o 2 TC 
1 = --l m trfr o V(p o = M,p = 0)]=a(N,N f )?-^ , (2) 

where a(3,2) = 1.40, o(2, 2) = 1.45 and a(2, 4) = 1.57. The same quantity 
was calculated in [I9| with the result a(3,3) = 1.43. A major problem with 



this calculation is that the gauge independence of 7 has not been established. 
With certain particular IR regularization schemes the result is gauge invariant 
but depends on how the IR regulator is removed and the external particle put 
on-shell [0, p. 

Moving p > massless fermion. 

There have been several papers concerning this case and most of them conclude 
that 7 depends on the magnetic screening E3, EJTL E3 . For abelian theories this 



is not an alternative so another resummation scheme was presented in [^5j (apart 
from the usual HTL resummation) which resulted in a non-exponential but finite 
decay of correlation functions. It is found that excitation with p ^> e 2 T (which 
is what we have) decay like 

e 2 T 

S{t) = e- %u * t exp[-v— tln(Mvt)\ , (3) 

47T 



where again v is the velocity and M. is the plasma mass. The expectation value 
of the damping rate is to leading order 

which has the same leading e 2 Tln - term as the heavy fermion case. The same 



leading damping rate has been obtained for the non-abelian plasma in |26, 



by assuming the existence of a magnetic mass of order g 2 T but the solution in 



|pq| is independent of such assumptions. Notice also that in [25] the constant in 
7 ~ e 2 T(ln - + const) was calculated. 

From the results above we conclude that, to the extent that it is calculable, the damp- 
ing rate of a light fermion due to scattering with gauge bosons is approximately 
7 ~ f° r both p = and p > 0. The correction factor ln(^p^) for moving 

fermions and the variation in the factor a(N, Np) is at most of order 2 and does not 
affect the conclusion significantly. 

For Higgsinos, the gauge boson contributions reflect weak and hypercharge symme- 
tries. For the contribution from the weak gauge bosons we obtain 7 oc 3(7|T/(167r) ~ 
0.0244T (g 2 = 0.640). The hypercharge contribution reads g\T /(16tt) ~ 0.00122T 
(gi = 0.247) which is negligible compared to the weak contribution. Thus we conclude 
that the gauge boson contribution to the damping rate of H ± is 

7gb ~ 0.025T (5) 

up to a factor 2 depending on non-perturbative effects through the magnetic mass and 
variations in the factor a(N, Nf). 

Let us now turn to scattering with scalars. There are a number of diagrams where 
the external fermion scatter with various scalars such as Higgs or squarks. They differ 
in structure depending on whether the fermion is Dirac or Majorana and whether the 
external Higgsino has a diagonal or off-diagonal mass term. However, the kinematics 
and coupling constants are all similar. Moreover, they do not suffer from the IR 
sensitivity of gauge boson scattering. Therefore, we shall only calculate explicitly 
the simplest diagram which is the stop-bottom loop with an external where the 
internal chiral fermion has no Majorana mass and the external particle has only a 
diagonal Dirac mass. Taking the notation from [E7[ the imaginary part is given by 



ImS P = - — / — — -sm20^-smh20 x _ P 

sin 2<pp J (27r) 4 2 2 

2n5((K-P) 2 -Ml)hl- l5 )2nA b (k ,k) , (6) 
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where 



A b (K) 



— (S b (k - ie, k) - S b (k + ie, k)) 
2m 

s(K)k j - r(K)k ■ 7 
s(K) 2 k 2 - r(K) 2 k 2 ' 



S b (K) 



r{K) 



s(K) 




(7) 



The notation is that Aib, M-h -^h± an d A* stands for the thermal bottom, stop, 
Higgsino mass and the vacuum Higgsino mass, respectively. Since s(K) and r(K) 
have imaginary parts for k% < k 2 there is a continuous part of the spectral function Ab 
apart from the 5-functions. We shall use the notation s(ko — ie, k) = Sr(K) + i si(K) 
and similarly for r(K). 

We want to take the expectation value of Im S in the external state of Higgsinos 
with thermally corrected masses. Since they have several branches and a continuous 
spectral weight below the light cone it is not possible to use a simple particle picture 
for the external states. One straightforward way to include collective effects is to 
multiply ImS with the spectral function for H^, i.e. Ag±(P), and take the trace and 
integrate over p Q . This gives the damping for a Higgsino averaged over all states with 
a given momentum. When the external states have several branches we have to be 
careful not to overcount the degrees of freedom, but the spectral function conveniently 
satisfies an exact sum rule which counts the number of degrees of freedom. Assuming 
that the external Higgsino mass is dominated by thermal effects we can use Eq. (|2|) 
for the damping rate: 



where we have divided by a factor 2 to average over the spins for incoming positive 
energy solutions. Finally we should do the K integral. For simplicity we choose p = 
so that there is no continuous contribution from the external states. Then the angular 
part is easy so we rewrite it as 





/ 



1 r°° 

d A K5{{K-P) 2 -M}) = — dk k{k )e{(k -po) 2 -Ml) . (9) 

07T° J-00 



Finally, we end up with 

1 r°° h 2 r 

p [ 3 >>0 H 



xsign(A;o - Po) 



+ 



X [ 0{kl ~ ^ ^k^^) S{K)k ^V {K ~ k ° ){k)) (10) 

+ 0(k 2 - k^ Slk °^ SR ~ ~ ^ R ~ r ^ k ^ ~~ 2s R k o[ s R s i k o ~ r R rik 2 } \ 



™ [(4 " ~ (rl - rf)Pp + 4[ SRS/ fc 2 - r^fc 2 ] 2 



where k — J (k — p ) 2 + A4^. This should be summed over the two solutions 



(±) 



'/i 2 + 4A^^ ± ± /i )/2. In first term of the curly bracket the 5-function is used for 
the ko integral so in practice one has to find all simultaneous solutions to 

Oo -Po) 2 - k 2 - M\ = , and D b (K) = , (11) 

and use the values of ko and k in the integrand. For the continuous part it is in 
principle necessary to integrate over all fco and for each value of ko and solve k from k = 



ko — Po) 2 — At the same time there are the factors 0((k o —p o ) 2 — M 2 ) 6{k 2 — kl) 
that restrict the integration domain. 

We have evaluated the damping rate for a Higgsino at rest numerically using ti~ = 
1 for the Higgsino Yukawa coupling (see Fig. ([I])). The scalar contribution to the 
thermalization rate (in the units of temperature) is presented for two different values 
of the effective thermal Higgsino mass: po = LOT (solid line) and po = 0.75T (dotted 
line), and for range of the right-handed stop mass form to T. It shows that for 
large enough stop mass (Ait > 0.5T for p = LOT and M-i > 0.25T for p = 0.75T) 
the main contribution to the damping rate comes from the (weak) gauge bosons. For 
the lower end of the stop mass range, however, the scalar contribution is considerably 
larger (7^ ~ 0.01T) but still smaller than 7 gb . 

The aim of this paper has been to establish a reliable value for the damping rate 
of Higgsinos in MSSM above the electroweak phase transition. Our claim is that the 
damping is typically dominated by scattering with gauge bosons and that the contri- 
bution from scattering with scalars is relatively small. We found that the gauge boson 
contribution to the damping is given by 7 g b ~ 0.025T. Scattering with stop/bottom 
gives a largest contribution for small stop mass (see Fig. (0)) j b a < 0.01T, but is 
much smaller for higher stop mass. Other scattering processes (involving gauginos 
and scalars Higgses) can be expected to give similar phase space integrals as the 
stop/bottom scattering but with gauge coupling constants instead of hi. We can there- 
fore conclude that they are relatively insignificant. In a similar way we can estimate 
the damping rate of H° and W. Actually, neutral Higgsino damping rate coincides 



Ib/t 




Figure 1: Scalar contribution to the Higgsino damping rate as a function of 
right-handed stop mass in units of T. The different curves are for po = LOT 
(solid line) and po = 0.75T (dotted line). 



with the charged Higgsino rates, because we are working in the unbroken phase. They 
also scatter with gauge bosons and we expect this to be the dominant process. In that 
case the damping rates would be 7^-0 — 7 9 f> — 0.025T and 7^ ~ g>|T /2n = 0.065T. 

We would like to reiterate why we only considered one single diagram including 
scattering with scalars. Even though all of them can be calculated with high accuracy 
using resummed propagators there is no reason why they should differ significantly 
numerically, apart from have different overall coupling constants. On the other hand, 
scattering of moving fermions with gauge bosons is highly non-trivial problem that can 
only be estimated using assumption about the value of the non-perturbative magnetic 
mass. Different estimates give qualitatively the same result but the accuracy can of 
course not be guaranteed. (The situation is different in abelian theories where there 
exists a resummation procedure controlling the IR sensitivity P5| .) It is, therefore, 
for the moment, no point in trying to achieve higher accuracy in the scalar scattering 
diagrams which anyway are sub dominant. The case is different for B° which is gauge 
singlet and thus does not have (one-loop) gauge boson contribution. However, the 
form of its scalar contribution differs from the scalar contribution of the Higgsinos 
because of the Majorana nature of B° |[28|| . 



Let us briefly discuss the implications of our findings. As we already mentioned, the 
precise knowledge of the thermalization rate of the supersymmetric particles is a key 
ingredient for the computation of the final baryon asymmetry. Sizeable decay rates of 
the particles propagating in the plasma destroy the quantum interference out of which 
the the CP-violating sources are built up and therefore reduce the baryon asymmetry. 
Small decay rates, on the other side, are relevant when the particles reflecting off the 
advancing bubble wall have comparable masses and resonance effects show up ||. In 
such a case, the thermalization rates provide the natural width of these resonances 
and as the present calculation demonstrates, in supersymmetric theories these depend 
in a complicated way on the particles involved and their plasma masses. Moreover, as 



shown in [g^, |30|, memory effects proportional to the damping rate of the Higgsinos 
play a crucial role in determining the lowest value of the CP- violating phases necessary 
to account for the observed baryon asymmetry. In particular, memory effects introduce 
an enhancement factor in the final baryon asymmetry of the order of 

Uhr) O 

with respect to the results found in ||. Here v w is the velocity of the bubble wall 
propagating in the plasma. Combining our findings with the analysis of Ref. |30 
indicate that the phase M of the parameter /i satisfy the inequality 

-2 l' V " 



|sin(0,)|>5xl(T^j. (13) 

These small values of the phase M are consistent with the constraints from the electric 
dipole moment of the neutron even if the squarks of the first and second generation 
have masses of the order of 300 GeV. 
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